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Abstract 

We propose a power filter Gp for linear reconstruction of the CMB signal 
from one-dimensional scans of observational maps. This Gp filter preserves 
the power spectrum of the CMB signal in contrast to the Wiener filter which 
diminishes the power spectrum of the reconstructed CMB signal. We demon- 
strate how peak statistics and a cluster analysis can be used to estimate the 
probability of the presence of a CMB signal in observational records. The ef- 
ficiency of the Gp filter is demonstrated on a toy model of an observational 
record consisting of a CMB signal and noise in the form of foreground point 
sources. 
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1 Introduction 



The exploration of anisotropy and polarization of the Cosmic Microwave Background 
(CMB) is fundamental in developing our knowledge about the Universe. During the 
last ten years several CMB measurement projects are being carried out ^. After the 
successful COBE mission the attention has been focused on obtaining higher angu- 
lar resolution of observational data in the range of the so-called Doppler peak and 
subsequent peaks. The level of the CMB anisotropy at multipoles / > 30 up to 
Z ~ 10^ — 2 • 10^ is a gold mine of cosmological information about the early Uni- 
verse and its most important parameters. However, statistical analysis of modern 
observational data is extremely complicated and time consuming. This complexity 
leads to development of different approximate methods of the data reduction such as 
Wiener-filtration radical compression method ^, likelihood method, band-power 
method and others The existence of a large number of processing methods of 
observational data for different types of experiments is easily understood. Firstly, 
there arc many difficulties and peculiarities of the measurements connected with the 
scan strategy, specific beam chopping, map making, etc. Secondly, the reconstruction 
of the parameters of the CMB signal is a sort of "inverse problem" . These difficulties 
are not specific for CMB measurements only. They are well known also in optical as- 
tronomy and in image reconstruction methods from satellite or airplane photography 
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Which one of the data reduction methods mentioned above is the most preferable 
for any type of experiments? The answer of this question depends on the our "intu- 
itive" hopes. We can use the following criterion which was suggested by Tegmark ^: 
"one method is better than another if it retains more of the cosmological information 
which operationally means that it will lead to smaller error bars on the parame- 
ter estimates". However, before CMB image reconstruction from the observational 
scans is done we have no information about the structure of the signal (systematic 
errors, Gaussianity, different types of foreground sources and noise, etc). So, roughly 
speaking, the Tegmark criterion plays an important role a posteriori, when the recon- 
struction of the CMB by different methods is performed. We have to remember that 
for different strategies of observations and for different concrete experiments different 
methods of data processing can be considered as the most preferable. In addition dif- 
ferent methods are preferable for different goals for using obtained data. In general, 
we may conclude that different methods complement each other. 

The purpose of this paper is to propose and investigate a new method with the 
following aim: how can we estimate the probability of the presence of the CMB signal 
in an observational map and determine its characteristics if we have a hypothesis 
about its statistical properties (for example that it is Gaussian) , and have a hypothesis 
about its power spectrum. 
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We suggest a new type of linear filtration of one- dimensional scans of CMB maps 
that is sensitive to non-Gaussian noise of any origin. This filtration does not change 
the primordial CMB power spectrum in contrast to the well known Wiener filter 
method. We use geometrical and topological characteristics such as the Minkowski 
functionals and statistics of peak distribution above and below some definite 
threshold ^2-i5_ show, that these characteristics are very sensitive criteria of 
Gaussianity of the restored signal. The filtration of the two-dimensional maps of the 
CMB with the filter which preserves the statistical properties of the underlying signal 
was described in 

The full advantage of the new method can be demonstrated even on very small 
data sets. We apply our method to a "toy" model of a mixture of CMB signal and 
foreground point sources to model balloon- and ground-based measurements at the 
frequency band u ~ 10 — 200GHz. The application of the method to real observational 
data will be considered in a separate paper. 

The rest of this paper is organized as follows. In Section 2 we describe the model of 
the real experiments and describe a new so called "power filter" for reconstructing the 
CMB signal. Section 3 is devoted to peak statistics and cluster analysis of the CMB 
signal in one-dimensional records. In Section 4 we present the new reconstruction 
technique of the CMB spectrum from records with point sources and demonstrate 
the strength of our method. Section 5 contains our conclusions. 

2 The model of real experiments and reconstruc- 
tion of the CMB signal 

The statistically isotropic distribution of the CMB temperature fluctuations on the 
sky is usually described by the spherical harmonics expansion Yim{q)'- 

oo I 

T{q) = ^Y1 aimBimYimiq), (1) 

(=1 m=—l 

where aim are zero-mean Gaussian deviates, {aima*,^,) = CiSu'Smm', with Ci given 
by the CMB power spectrum, Bi^, the antenna beam, and q{6, ip) the unit vector 
defining the position of the observational direction on the sky. Eq. (|I]) can be written 
for all frequency channels z/i...z/m, where fh is the number of bands. After the COBE 
mission, different CMB experiments such as Saskatoon, QMAP, CAT and others 
have observed rather small parts of the sky with higher angular resolution. They 
estimated statistical characteristics of the CMB signal from the maps using some 
specific pixelization strategy. Two future space missions, MAP and PLANCK, will 
cover a significantly greater part of the sky and will provide data with higher precision. 
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The traditional view is that most information of a Gaussian random field is con- 
tained in the multipole power spectrum Ci. As a result a great deal of efforts has 
gone into obtaining the best estimate of the power spectrum Ci from the experimental 
data. The power spectrum Ci describes properties of a two-dimensional map. How- 
ever, there exist many balloon-borne experiments and ground-based projects which 
performed observations of the CMB anisotropy along one- dimensional sky patches. 
Besides that the future MAP and PLANCK missions will collect data from a large 
number of intersecting circles which will then be reconstructed into a two-dimensional 
sky map. Briefly speaking, because of the simplicity of data analysis in the one- 
dimensional case, the collection of CMB observational data in this format is an in- 
teresting option for modern anisotropy experiments. 

Usually (see ref. ^ for a recent review) one has a time-ordered set of data and 
can, knowing the adopted scan strategy, obtain a map with pixel values denoting the 
sum of a pure CMB signal x and different types n of noise, including foregrounds, 
atmospheric emission, radiometer noise and so on. 

In the following we consider final maps consisting of a pure CMB signal x and 
noise n = and use therefore for a one-dimensional single record 

Y(t) = x(t) + n(t) (2) 

where t is an independent variable along the record. 

For this simple model of the measured signal Eq.(^ one can approximately recon- 
struct the pure CMB signal using a linear reconstruction technique (see ^ for a general 
description of the method). Note that in Eq.(0) we do not assume the random fields 
x and n to be Gaussian. 

The general task of linear methods of reconstruction of a signal is to obtain an 
approximate estimation x'(t) of a signal x(t) from the measured signal Y(t) using 
linear operations. Following ^ the estimation x'(t) of the CMB signal x(t) from Eq.(|^) 
can be written as 

x'{t) = J G{t-t')Y{t')dt' (3) 

where G{t — t') is the filter of the reconstruction. 

Before describing the new method which we propose for the reconstruction of 
the CMB signal, let us briefly discuss the standard linear filtering methods which 
were reviewed by Tegmark ^. Firstly, following ^, we point out that linear methods 
of filtration of maps do not destroy information about phases which is contained in 
the initial maps (or one-dimensional records). Secondly, the non-linear extraction of 
the signal, such as band-power method needs some special assumptions about 
the statistical properties of CMB signal and noise (for example Gaussian statistics 
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of noise). Besides that, all non-linear methods destroy information contained in the 
maps ^. 

Let us return to Eq. (j^) and Eq. (^) and consider one example of CMB signal recon- 
struction which uses a simple Wiener filtering method ^. Following ^, the estimated 
signal from Eq.(|]) is 

x"{t) = Jw{t- t')Y{t)dt, (4) 

where W{t — t') is the so called Wiener filter which minimizes the reconstruction error 
= (x" — x)^. Tegmark and Efstathiou ^ and Tegmark ^ have shown that the W 
filter has a particularly simple form in Fourier space under assumptions in Eq.(0). 
For this case W has the following form: 

,Y(JA _ PcMBjk) _ f Pnoise{k) \ ^ , . 

where Ptot{k) = PcMsik) + Pnoiseik) is the power spectrum of the measured signal 
Y(t), PcMB{k) the CMB power spectrum folded with the antenna beam and Pnoiseiki) 
the noise power spectrum. 

Let us return to Eq.(|). It is straightforward to prove that the power spectrum 
of the reconstructed signal Px"{k) is equal to 

Px"{k) = W'^{k)Ptot{k) = PcMB {k)W{k). (6) 

As one can see from Eq.(|^), W < 1 which, together with Eq.(|]), shows that the 
Wiener filter decreases the CMB power spectrum of the reconstructed signal. 

For a linear reconstruction of the CMB signal in one-dimensional scans we propose 
to use the so called "power filter" Gp ^. As we told in Section 1 the corresponding 
consideration for two-dimensional CMB maps was done in ^. This filter minimizes 
the difference between the power spectrum Px'{k) and the CMB power spectrum. It 
is clear that this condition corresponds to 

Px'{k)=PcMB. (7) 

We will demonstrate that for our case Gp = W^^'^. For the proof let us rewrite Eq.(|^) 
for our filter Gp-. 

x'(t) = j Gp{t - t')Y{t')dt'. (8) 

Now note that from Eq.(§) one can conclude that 

Px'ik) = Gl{k)Ptotik). (9) 
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Using Eq.(ll), Eq.(|^) and Eq.(P) we conclude that 



GIW = = (10) 

Ptot{k) 

which proves our statement: Gp = W^^"^ . 

3 Peak statistics and cluster analysis of the recon- 
structed CMB signal 

Using power fihration by a Gp filter we transform the initial observational map to a 
reconstructed record {x'(t)}. Like any other linear filter the Gp filter transforms the 
power spectrum of the initial signal preserving phases of the signal. This property 
of the linear filtration is extremely important and we will consider it in more details. 
Let us consider a one- dimensional record Y{t) and its Fourier expansion: 



°° f2Tmt\ , f2Tint\ 

an COS ( I + On siu 



n=l 



(11) 



where T is the total length of a record. After a power filtration we get 



°° ' ^ f27mt\ , ^ . /27mt\ 
a„G„ cos —~— + bnGn sm 



n=l 



(12) 



The initial power spectrum of the record (Eq.(ll)) is given by 



Pt^ = lial + bl) (13) 
and of the reconstructed power spectrum (Eq.(12)) by 

Pt''^ = lGl{al + b'J. (14) 
However our Gp filter does not destroy the spectrum of phases since 

= arctan | — ) = (15) 

So, roughly speaking, we transform the power spectrum of the initial signal by Gp 
filter but preserve the spectrum of phases. Thus Eq.(|l^) demonstrates that such a 
filtration preserves a Gaussian or non-Gaussian property of the modes which are not 
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suppressed by the filter. This property is vahd for all linear filters. The specific prop- 
erty of the Gp filtration is that in addition to preserve the phases it also preserves the 
spectrum of the CMB signal. Gp is a unique filter with this property. Cluster analysis 
and peak statistics ^^"^^ are very sensitive to both phases and the power spectrum. 
Thus, if we believe that after the Gp filtration the reconstructed signal really repre- 
sents the approximate distribution of AT of CMB signal, it should have the same 
phase characteristics and the same peak statistics and properties of clusterisation 
which are predicted for the initial CMB signal (under an accepted hypothesis about 
the cosmological model). Confrontation of peak statistics and clusterisation of the 
reconstructed signal with the theoretical predictions is a good test of the "quafity" 
of the reconstruction. 

We will investigate how this method works in the case of one-dimensional records. 
In the previous section we have shown that the reconstructed record x'(t) has the 
power spectrum Pcmb- We would like to point out that all non-Gaussian features 
of the initial signal are transformed to the reconstructed record x'(t) because of the 
finear character of filtration. We also note that in the case of application of the Wiener 
filter, peak statistics of the reconstructed signal is different from the hypothetical 
CMB signal because of the change of the power spectrum. 

Assuming the primordial cosmological signal to be Gaussian, we perform the test 
of the filtered signal with regard to Gaussianity. In fact, we calculate the probability 
that the filtered signal has a Gaussian nature. 

The Gp filtering obviously makes the power spectrum of the observed signal equal 
to the expected CMB spectrum. The rest of the information to be investigated is 
about the phases only. In order to perform this test, we use a geometrical approach. 



We normalize the x'(i) record to its variance: v{t) = AT/y'(AT^) and consider 
the following characteristics of the v{t) statistics: 

A) Normalized distribution functions of maxima and minima above the threshold 

Jmax\yt) — ^ Jmin\yt) — ^ , Ut)j 

where Nmax{^t) and Nmin{^t) are the number of maxima and minima, respectively, 
and Nmax = Nmax{n ^ "Oo), Nmin = Nmin{'t^t ^ "Oo) are the total number of 
maxima and minima. 

B) The global Minkowski functionals Mi and M2 for a one-dimensional record, 
which characterize the cumulative distribution function (Mi) and the Euler charac- 
teristic (genus) of the random signal (M2). The last quantity was first introduced in 
cosmology by Doroshkevich For the N dimensional case there are N-l-1 Minkowski 
functionals. 
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C) Finally we consider the mean length of the clusters above and below a constant 
threshold Ut. A one-dimensional cluster of maxima is the continuous part of the curve 
with i/(t) > z/j inside the interval t G (ti, ^2)- Here ti and ^2 are roots of the equation 
^{^1,2) = and t2 > ti. The length k{yt) of the cluster is defined as the number of 
peaks with height Vpeak > in the cluster. If the value of z/^ is high (z/^ ^ 1) then 
only high maxima (separated from each other by regions with v{t) < ut) are present 
above the threshold and the typical length of a cluster is {k{vt)) = 1- The reduction 
of the threshold level z/f down to z^t — > or z/^ < leads to the appearance of big 
clusters when maxima of smaller clusters begin to connect together and generate a 
new cluster. This process is characterized by the number of clusters Nkiyt) of the 
length k and the total number of clusters N{vt) which are present in the anisotropy 
records for an appropriate threshold z/^: 

N{v,) = T.t=,N,{vt). (17) 

Thus for different statistics of initial records the mean length of a cluster at threshold 
level Vt is 



It is necessary to note, that only the first characteristic considered under A and the 
first Minkowski functional Mi are independent. The rest of them can be represented in 
terms of the distribution of peaks (A). Thus, we use them only for visual clarification 
and representation of our results. 

The defining Eq.(^) - Eq.([TE|) do not depend on the nature of the initial signal 
or its composition. 

Below we use the standard definition of the spectral parameters: 

al = je^PcMB{k)dk, n = 0,1,2 

a-oo-2 

where PcMsik) is the one-dimensional CMB power spectrum. 
The characteristics A, B and C have a rather simple form in case of a random 
Gaussian field 10-12,14. 

The number of maxima and minima above a threshold z/j are: 



N = N ■ = — ^ 

' max ^ ' mm 2-k ai 



^2(1-7^)/ ^ V VV2{T^ 



N^U^t) = (1 - $ f^^^) - 7e-^ (l-^( ^^^'I'l'l (2°) 
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The Global Minkowski functionals for a one-dimensional Gaussian field are given by: 

(21) 



2 
It 



where $(x) ^^J^e-''dt. 

As it has been already mentioned above, the last characteristic (M2) can be written 
in terms of the peak distributions (A) 

It is not difficult to demonstrate that for a Gaussian distribution all these character- 
istics are functions of only two variables: threshold Vt and spectral parameter 7. This 
fact makes our analysis transparent. 

To calculate the probability that the restored (filtered) signal is Gaussian, we 
use only the characteristic as mentioned under A. The distributions of maxima and 
minima, fmaxi^til) a-nd fmin{j^t,l), above a threshold Ut are monotonic functions of 
i/f with 

fmax{ OO) T) — fmin{ OO) T) — ^ 
fmax(+00, 7) = fmin(+00, 7) = 0. 

For a limited number of peaks these distributions look like a step functions with num- 
ber of steps equal to the total number of maxima and minima. If the observed signal 
has a cosmological nature, then this function should be consistent with the smooth 
analytical curve (see Eq. (16) and (20)) for 7 of the CMB spectrum. In case of real 
observations in a limited record we deal with bad statistics (small number of peaks). 
For this kind of statistics we have to perform an appropriate test, e.g. a Kolmogorov- 
Smirnov test Before applying such a test and to demonstrate the advantage of our 
method, wc should do the following: 

1. As representations of the spectrum of the filtered signal (being identical to the 
spectrum of the hypothetical CMB signal) we simulate a large number rir of different 
reahzations of a random Gaussian field in the considered record using a FFT algo- 
rithm. 

2. For each realization r=l,..,nr. we calculate 7^ using the definition (19) as well as 

and using Eq.(16). Furthermore we calculate the quantities: 

Imax,real{r,'yr) = S-5 fmax{^Ulr)dUt 

Imax,anal{j'i'^r) /— 5 fmax{yti'^r)d'Ut (24) 
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and 



^max -^maXjreal -^max,anal (0K\ 

A . — T . . _ r . . y"^^) 



where fmaxi (z/t, 7^) have the analytical form (20) and f^axi^u 7r), fLni^t, % 

are distributions of maxima and minima in a given realization. In (24) we performed 
integrations with the limits (-5, +5) where ]raax and /^m in practice have reached their 
asymptotic values given by Eq.(23). Obviously, 7 will differ slightly from one realiza- 
tion to another because of the cosmic variance. This difference becomes smaller for 
longer records. Now we derive the formal probability that the distribution of peaks 
corresponds to a Gaussian distribution. There are different methods to do this. One 
of them is the following: 

3. We count the fraction of realizations, where |ATOaa;(r)| and |Ami„(r)| are greater 
than some positive value z and calculate P{z) 

p (^\ _ N{\Amax{r)\>z) 
-' max \^ ) 



r 



N{\/^J'^(r)\>z) (26) 

min\^J „ 



Here P{z) represents a numerical estimate of how well the distribution of maxima 
and minima corresponds to a Gaussian signal. 

4. For a given filtered record we can now estimate to what extent the signal has a 
Gaussian nature. We derive A^nax ^-nd A^j„ to determine the probability P{z). It 
is clear that if P{z) is larger than ^ 3% and therefore z is within the 3a level of 
the distribution, there is a good chance that we observe a Gaussian signal. Alterna- 
tively one can calculate the differential functions n(A) of the distributions of A^nax 
and Ajnin and consider the actual value n(A) for the filtered signal as a measure of 
closeness to a Gaussian signal. If this n(A) is greater than a few percent the filtered 
signal is Gaussian within the " 3a level" using the analogy with a standard Gaussian 
distribution. 

One could also use a Kolmogorov-Smirnov test, some of its generalizations or other 
methods to be discussed in a separate paper. 

In the next section we consider the application of this technique to a concrete toy 
model of an observed record. 



4 Application of the power filtration to a toy model 
of an observational record 

Now we demonstrate the application of the Gp filtration and estimate to what extent 
we can recover the Gaussian nature of the signal. Furthermore we compare with the 
result from a Wiener filtration. 
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We consider a CMB spectrum of a standard Cold Dark Matter cosmological model 
with parameters: /iioo = 0.75, Q = 1, Qb = 0.0125, Qa = 0. For this spectrum we 
simulate = 1000 different realizations of the random Gaussian field (using a FFT 
algorithm) for a record with length 360° and with FWHM of the antenna beam equal 
to 0.5°. Using the definitions in the previous section we calculate the values of 7^, 
Imax, Imin ^oi cach realization. These values are slightly different from one record to 
another because of the cosmic variance. In Fig. 1 we show the dependence of I max 
and Imin on 7^ for the realizations. In the same Fig. 1 we also show the corresponding 
Gaussian dependences. 

Now we consider one specific realization (in analogy to a realization on the sky) 
which is marked by a star on Fig. 1 (7= 0.3738, Imax= 2.751 and Imin= 2.271). Next 
we compute for this realization f{i't) (see Fig. 2a), Mi(z/t) (see Fig. 2b), rate of 
clusterisation {Nmax/ {Nmax - Nmin) and Nmin/ {Nmin - Nmax) as a function of ft (see 
Fig. 2c), and M2{i't) (see Fig. 2d). On Fig. 2(a-d) we also show the corresponding 
analytical dependence (see Eq.(PD|) and Eq.(pTD). The graphs for the realization on 
Fig. 2(a-d) are only slightly different from the corresponding analytical curves given 
by dashed lines. 

Then we add to the CMB signal a signal from 100 identical foreground point-like 
sources with a random distribution over the record and with an amplitude such that 
the rms of the maps increases by a factor \/2 relative to the pure CMB signal. We 
consider the signal from the point-like sources as "noise". The signal/noise ratio in 
this example is thus equal to 1. Of course this is not a representation of a realistic 
noise but only used to demonstrate the efficiency of the power filter Gp. Discussion 
of more realistic models will be given in a separate paper. The power spectrum of 
the CMB signal plus noise is computed for the same antenna beam with FWHM 0.5°. 
For this realization we determine ryC^B+noise Imax and Imin- The corresponding 
points are marked by crosses on Fig. 1 for the parameters 7= 0.4833, Imax= 2.871 
and Imin= 2.248. 

It is seen that the crosses are far away from the clouds of realizations of the CMB 
and far from the Gaussian points (marked by dots with the same ^CMB+noise pjg_ 
1). For this realization we calculated fmax, fmin, Mi, the rate of clusterisation, and 
M2 and the corresponding analytical Gaussian functions for the given ^CMB+noise^ 
The corresponding plots are presented on Fig. 3. One can see that the distributions 
of these values on Fig. 3(a-d) differ significantly from the corresponding Gaussian 
functions. This means that the CMB+noise signal is far from being Gaussian. The 
same is seen from Fig.l. Now we perform filtration of this CMB+noise signal using 
the power filter Gp. As a result we expect to obtain a reconstructed signal with 
the power spectrum nearly identical to the CMB power spectrum of the considered 
realization. The corresponding points for I^ax and Imin are indicated in Fig. 1 by 
circles with the parameters 7= 0.3930, Imax= 2.782 and Imin= 2.272. 
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It is seen that the circles in fact are within the corresponding clouds of numerical 
realizations of the CMB signal and close to the corresponding Gaussian points with 

the same 7. 

In Fig. 4 we have plotted fmax, fmin, Mi, the rate of clusterisation, and M2 for 
the reconstructed signal as well as the corresponding Gaussian analytical functions. 
One sees that the reconstructed realization is considerably closer to Gaussianity than 
the CMB+noise signal, demonstrating how effective the power filter works. It is 
particularly important that the filtered signal has practically the same 7 as the input 
CMB signal. Furthermore we may get a numerical estimate of how good the power 
filtration really is by using the method described in Section 3. The values of n(A) 
for the filtered signal are n{Amax)— 0.12 and n(A^j„)= 0.51. 

For comparison we have performed a Wiener filtration of the same realization of 
CMB+noise signal and computed 7*^, and The values (7= 0.6232, Imax= 
2.930 and Imin= 2.067) are plotted as triangles on Fig. 1 together with the Gaussian 
values for the same 7^ as dots. 

One can see that the Wiener filter displaces the corresponding points far away 
from the clouds of realizations of the CMB signal and that the power spectrum of the 
reconstructed signal is very far to be spectrum of the CMB signal. 

5 Conclusions 

In this paper we propose to use a power filter ^ for linear reconstruction of the CMB 
signal in one-dimensional scans. The specific property of the filter is that it 
preserves the power spectrum of the CMB signal. Under conditions described in 
Section 2 we get Gp = W^^"^ where W is the Wiener filter. Unlike the power filter 
Gp the Wiener filter W diminishes the power spectrum of the reconstructed signal 
compared to the CMB power spectrum. In section 4 we have demonstrated that 
peak statistics and cluster analysis can be used to estimate the efficiency of the 
power filtration and to estimate the probability that a CMB signal is present in an 
observational record. In Section 4 we also demonstrated how the power filter works 
using a toy model of an observational record consisting of the CMB signal and noise 
in the form of point like sources. We note that in the example of Section 4 we used 
a signal/noise ratio equal to unity as described in section 4, but tests with different 
toy models have demonstrated that this method works well also in the case with this 
ratio significantly less than unity. Discussions of more realistic models of the noise 
as well as application of the described method to real observations will be done in 
subsequent papers. 
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Figure 1: Dependences Imax (upper part on vertical axis) and Imin (lower part 
on vertical axis) of 7 (horizontal axis). Thick dots are dependences for Gaussian 
analytical expressions (20), (24). Stars correspond to the CMB realization with 
{1 J max, Imin) = (0.3738, 2.751, 2.271). Crosses correspond to the CMB+noise re- 
alization with {'jJmaxJmin) = (0.4833, 2.871, 2.248). Circles correspond to the re- 
constructed signal using the Gp filter, parameters are {j, Imax, Imin) = (0.3930, 2.782, 
2.272). Triangles correspond to the reconstructed signal with the Wiener filter; pa- 
rameters are {■y, I max, Imin) = (0.6232, 2.930, 2.067). 
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Figure 2: Dependences a) fmax, fniiiu b) Mi, c) the rate of clusterisation (k) and 
d) M2 as function of Uf for a CMB realization marked on Fig. 1 by stars. Dashed 
lines correspond to analytical Gaussian expressions (see text), fmax and fmin are 
normalized by N^ax + Nmin- 
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